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ABSTRACT 

We present a new technique for determining an upper limit for the mass of the black 
hole in active galactic nuclei showing warm absorption features. The method relies on 
the balance of radiative and gravitational forces acting on outflowing warm absorber 
clouds. It has been applied to 6 objects: five Seyfert 1 galaxies: IC 4329a, MCG— 6-30- 
15, NGC 3516, NGC 4051 and NGC 5548; and one radio-quiet quasar: MR 2251-178. 
We discuss our result in comparison with other methods. The procedure could also 
be applied to any other radiatively driven optically thin outflow in which the spectral 
band covering the major absorption is directly observed. 

Key words: galaxies: active — galaxies: Seyfert — X-rays: galaxies. 
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massive black liulu (OMDII) (see pjy udmi-Bell 100 



and |Re es 1984). For many years, theorists and observers 



have attempted to prove the existence of SMBHs, inferring 
their presence by studying how they affect t he su rrounding 



environmen t (Kormendy & Richstone 1995 and Richstone 
| et al. 1998 review the SMBH search). A variety of tech- 



niques and arguments has been de veloped to estimate th e 
mass of the central engine in AGN ( Ho & Kormendy 2000 ) . 



This has included both theoretical and observational meth- 
ods such as: fitting the spectra of accretion discs; high ef- 
ficiency and Eddington limit arguments; stellar dynamical 
searches; kinematics of radio masers; X-ray variability and 
virial masses from optical variability, which includes rever- 
beration mapping and photoionization model techniques. 

In this paper we propose the use of warm absorbers 
(WAs) to determine an upper limit to the BH mass. WAs are 
significant quantities of partially-ionized, optically-thin ma- 
terial found mainly in the vicinity of Type 1 Seyfert galax- 
ies. This gas is most readily observable when it lies along 
line-of-sight to the central X-ray source: the observed X- 
ray sp ectrum then d isplays well-defined K-shell absorption 
edges. Halpern (1984) was the first to postulate the presence 
of such material based on the unusual X-ray spectrum of 
the QSO MR 2251-178 from the Einstein observatory. The 
clearest demonstration of the presence of warm absorbing 
material comes from spectroscopic observations of the OVII 
and OVIII K-shell absorption edges at 0.74 and 0.87 keV , 
respectively. These edges were first found by the ROSAT 



et al. 1997 



PSPC ( |Nandra fc Pounds 1992j ). The spectral capabi lities of 
ASCA allowed for detailed studies of such edges (e. g. Fabian 
et al. 199^, | Otani et al. 1996j prandt et al. 1996; , |Mathur 
eynolds 1907. |Lcighly et al. 1997 and George 



et al. 199j ). WAs are common, being detect ed in half of the 
sources " peynolds 1997] , |George et al. 1998[). Several differ- 



ent models have been suggested concerning the nature and 
location of the W A. Among others, these incl ude the very de - 



tailed studies by Krolik & Kriss (1995) and Netzer (1996) 



The reason for the large variety of models is that not all 
physical properties of the WA can be directly determined 
from X-ray spectral fits, but only certain combinations of 
them. 



In addition to the WAs, associated UV and optical 
absorpti on features are also t ypically found in AGN spec- 
tra (see |Crenshaw et al. 1999 and references therein for an 
overview). A common property shared by almost all associ- 
ated absorption features is that any velocity shift is always to 
the blue of the systemic velocity, implying that the absorbing 
material is outflowing, prenshaw et al. (1999 ) found a one- 
to-one correspondence between those objects showing UV 
absorption and those presenting X-ray WAs. This implies 
the two phenomena are likely to be physically related. How- 
ever, since WAs generally reveal higher ionization states, 
and higher column densities, they can not be the same in 
general. A very plausible connection is a dynamical one in 
which both UV absorbers and X-ray WAs are outflowing. 
The blue-shifts in X-ray absorption line s have recently been 
revealed by Chandra observations (e .g. Kaastra et al. 200C , 
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Kaspi et al. 2000, dollinge et al. 2001), presenting typical ve- 
locities of several hundred to several thousand km s _1 . One 
unresolved issue is whether the outflowing gas is smoothly 
distributed (i.e. winds) or resides in discrete clouds. Krolik & 
Kriss (1995|) show that the warm absorption seen in many 
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X-ray spectra of AGN may be produced by X-ray heated 



winds. Reynolds & Fabian (1995) also sketch some ideas re- 
lating the WA to a radiatively driven outflow . In contrast to 
the wind model, Shelouche & Netzer (2001) show that ac- 



celerations to velocities of a few hundred to a few thousand 
km s _1 is a natural consequence of the AGN environment 
when radiation pressure acts on external pressure-confined 
highly ionized clouds (responsible for the WA features). For 
the present study, we will adopt the latter model in which 
WA clouds are radiatively accelerated by the central ionizing 
continuum. 

This work describes a new method for obtaining an up- 
per limit on an AGN BH mass, by means of a comparison 
between the radiative and gravitational forces acting on its 
WA clouds. Section 2 presents the 30 observations of the 6 
AGNs sample to which our technique has been applied. The 
description of our method, a spectral analysis of the A SCA 
data and the results and comparison with other BH mass 
estimates by different techniques are addressed in Section 3. 
Section 4 and 5 analyse the influence of the underlying con- 
tinuum on our estimate and further terms in the equation 
of motion, respectively. Finally, in Section 6 our conclusions 
are summarized. Throughout this paper, unless otherwise 
stated, errors on physical quantities are quoted at the lcr 
level for one interesting parameter. 



2 THE DATA 

We have studied a sample of six AGNs consisting of five 
Seyfert 1 galaxies (IC 4329a, MCG-6-30-15, NGC 3516, 
NGC 4051 and NGC 5548) and one radio-quiet quasar (MR 
2251-178), which cover a range of almost 4 orders of mag- 
nitude in X-ray luminosity. These objects are a subsam- 
ple of those presenting warm absorption features in the 24 
Type 1 AGN sample of Reynolds (1997). With the excep- 
tion of MR 22 51-178, the other five Seyfert galaxies are also 
present in the George et al. (1998) survey. The WAs o f all 



these source s have a very rich observational history (Ko- 



mossa 1999) and their host AGNs are amongst the bright- 
est ones in the nearby Universe (z > 0.07). The spectral 
files have been downloaded from the Tartarus Database 



( xttp : //tartarus .gsfc .nasa.gov/) and the 30 observa- 
tions analysed are listed in Table 1. 



3 UPPER BLACK HOLE MASS ESTIMATE 
3.1 Description of the method 

The following assumptions have been made: 

• the WA, responsible for the absorption features around 
1 keV , is present in the form of clouds. 

• these clouds are pushed radiatively outwards. 

• the clouds motion is entirely radial. 



The radial equation of motion for these WA clouds is 
AL ahs GMbhtti 



■ p v A = m wa a 



(1) 



4irr 2 c r 2 
where: A is the cross-section presented by the WA clouds; 
Labs is the luminosity absorbed by the WA cloud; r is the 
distance from the ionizing source to the WA cloud; c is the 



Source 


Type 


Redshift 
z 


Galactic Njj 
(10 20 cm,- 2 ) 


Date of 
Observation 


IC 4329a 


Syl 


0.016 


4.55 a 


93/08/15 
97/08/07 
97/08/10 
97/08/12 
97/08/15 


MCG-6-30-15 


Syl 


0.008 


4.06 a 


93/07/09 
93/07/31 
94/07/23 
97/08/03 
97/08/07 


MR 2251-178 


RQQ 


0.068 


2.8 


93/11/16 
93/12/07 
93/12/14 
93/12/19 
93/12/24 


NGC3516 


Syl 


0.009 


3.4 


94/04/02 
95/03/11 
95/03/12 


NGC4051 


Syl 


0.002 


1.31" 


93/04/25 
94/06/07 


NGC5548 


Syl 


0.017 


1.7 


93/07/27 



96/06/27 
96/06/29 
96/07/01 
96/07/03 
96/07/04 
98/06/15 
98/06/20 
98/07/07 
99/01/19 



Table 1. The AGN sample. Column 1 gives the common source 
name. Column 2 indicates the type of nuclear activity (Syl = 
Seyfert 1 galaxy; RQQ = radio-quiet quasar). Column 3 give 
the redshift of the sources (Veron-Cetty & Veron 1993). Col- 
umn 4 gives the Galactic HI column density towards the source 
as determined by 21-cm measurements ["a" indicates an accu- 
rate value obtained from Elvis, Wilkes & Lockman (1989); other- 
wise quoted values are interpolations from the measurements of 
Stark et al. (1992)]. Column 5 gives the year /month/day of the 
ASCA observation used (see the Tartarus Database for a detailed 
description of each observation). 



velocity of light; G is the gravitational constant; Mbh is the 
black hole mass; m wa is the mass of a WA cloud; p is the 
density of the medium through which the WA is moving; v is 
the velocity of the WA clouds with respect to that medium 
and a is the acceleration experienced by the WA cloud. 

The first term on the left hand side is the outward force 
due to radiation pressure, the second term is the restoring 
force due to gravity, and the third term describes any drag 
force (i.e. ram pressure) experienced by the WA when mov- 
ing through an intercloud medium. If the drag force and the 
acceleration term are negligible, a comparison between the 
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radiation and gravitation terms leads to an expression for 
the BH mass: 



47rcGm p jV wa 



(2) 



where m p is the proton rest-mass and iV wa is the column 
density of the WA clouds. If these terms are non negligible 
and positive, then the above expression becomes an upper 
limit for Mbh- We shall treat it as such in the rest of the 
paper and discuss further the drag and acceleration terms 
in Section 5. 

The only two unknowns in Eqn. 2 are L a t, s and iV wa . 
iV wa is calculated as a parameter in our spectral fits of a sin- 
gle ionization parameter mod el presented in the fo llowing 

Subsectj£lB As already seen in Reynolds (1QQ7 ) and ripriT-gr 



respectively. These two edges provide, to a good approxima- 
tion, a description of the effects of the warm absorber over 
the ASCA band. 

• intrinsic absorption by a column density Vh of neutral 
matter in the rest frame of the source. This absorbing ma- 
terial is assumed to have cosmic abundances. 

• Galactic absorption (at 2 = 0) fixed at the level deter- 
mined by the HI 21-cm measurements reported in Table 1, 
column 4. 

Free parameters in the fit were V, tot, ros and JVh- The 
normalizations of the model for each of the four instruments 
were also left as free parameters. This allows for the known 
~ 10-20 per cent discrepancie s between the normalizations 
of different instruments (see [l ttp : //heasarc .gsf c .nasa. 



et al. ( Jflfllj). the data can often be well fit bv such a model gov/docs/asca/watchout .html| for a discussion of the de- 



which allows for accurate estimates of the ionization pa- 
rameter and column densities. The calculation of L a bs has 
been performed using the fact that the photoionized plasma 
causing the edges is optically thin along the line of sight. 
Therefore, the expression for the radiation pressure, P ra( j, 
can be written as 



Prad — 



4 7rr 2 c 47rcr 2 



(L(u) - L ohs (v))dv 



(3) 



where: L(y) is the luminosity per unit frequency from the 
ionizing source (which we assume to be a power law of pho- 
ton index F) and L b s (y) is the luminosity per unit frequency 
observed after L(y) has passed through the WA cloud. The 
photoionization models can predict both the intensity of the 
incident and transmitted continuum allowing the calculation 
of the above integral (see Section 3.2.3). We have chosen the 
limits of our integration to be: v m i„ = 1 eVand i/ m „ = 20 
keV. 



crease in efficiency of SISO and SIS1 below 1 keV since ap- 
proximately late 1994). Table 2 presents a summary of the 
X 2 fitting. 

As noted in Reynolds (1997), fixing the energies of the 
two absorption edges at the physical rest-frame energies of 
the K-shell absorption edges of OVII and OVIII, is the most 
robust to the various unmodeled spectral complexities such 
as recombination line/continuum emission, resonance ab- 
sorption lines and other absorption edges. These complex- 
ities would require many additional degrees of freedom to 
model (such as the covering fraction of the warm material, 
its density, the chemical abundances and the velocity struc- 
ture) and would lead to an over modeling of the data for 
all but the very highest quality ASCAdata,. In fitting a sim- 
ple two-edge model, these complexities could lead to false 
shifts in the threshold edge of the edges (e.g. the existence 
of the Ne IX edge may falsely imply a small blueshift of the 
OVIII edge when a two-edge model is fitted with the edge 
threshold energies as free parameters). 



3.2 Spectral analysis 

The data have been fitted using two different models. The 
aim for this multiple analysis has been firstly, to determine 
phenomenologically the underlying continuum by mimick- 
ing the WA features by the two photoionization absorption 
edges at the OVII and OVIII rest-frame energies (Sections 
3.2.1 and 3.2.2); and then to calculate the radiation pres- 
sure acting on the WA (Section 3.2.3). Data from all four 
ASCA instruments were fitted simultaneously. SIS data be- 
tween 0.5-4.0 keVand GIS data between 0.8-4.0 keVwere 
used. The spectra were then fitted using the standard \ 2 
minimization technique implem ented in versio n 11.0.1 of the 
XSPEC spectral fitting package (Arnaud 1996). 



3.2.2 Soft excess 

A detailed study of the soft excess was not performed since 
soft excesses are expected to be most noticeable below ~ 0.6 
keV . This is the energy below which significant SIS calibra- 
tion uncertainties exist. However, for completeness and fol- 
lowing the analysis in Reynolds (1997), a black-body compo- 
nent (subject to neutral and ionized absorption) was added 
to the above model and the data were re-fitted. To avoid 
being severely biased by calibration uncertainties below 0.5 
keV , all data below this energy were discarded. The results 
found are presented in Table 3, where a significant improve- 
ment of the fit is obtained for MCG-6-30-15, NGC 3516, 
NGC 4051, and NGC 5548. 



3.2.1 Phenomenological model 

In order to determine the underlying continuum hidden by 
various spectral features displayed by the data, the spectrum 
of each object was fitted by a spectral model consisting of 
the following components: 

• a power-law representing the primary continuum (pho- 
ton index Y) 

• two absorption edges with rest-frame threshold energies 
of 0.74 keVand 0.87 keV representing OVII and OVIII K- 
shell absorption (with maximum optical depths ro7 and ros) 



3.2.3 Photoionization models 

It will be assumed that the physics of the WA is dominated 
by photoionization. This assumption results from the very 
high and variable ionization state. Quantitative evidence 
for this comes from the observed anti-correlation between 
the OVIII absorption edge depth and the primary ionizing 
flux observed during t he long observation of MCG— 6-30- 
15 jOtani et al. 1996| ). Also NGC 4051 has shown vari- 
ability of the OVII edge (Guainazzi et al. 1996). A grid 
of self-consistent models has been produced using Ferland's 
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Source 


Date of 
Observation 


Photon index 

r 


TOT 


T 08 


Intrinsic Nh 
(10 22 cm" 2 ) 


X 2 /dof 


Red. x 2 


IC4329a 


93/08/15 
97/08/07 
97/08/10 
97/08/12 
97/08/15 


1 90+ 01 
1,w -0.01 
I yn+O.Ol 
1 - ' a -0.02 
1 ni+0.02 
1 - al -0.02 
1 Ql+0.02 
^ -0.02 


67+°' 0e 
u -°' -0.06 

11+ 007 

u ' x -0.07 
45+O.O6 
u -^°-0.08 

n 90+0. 08 

0QQ + 0.10 

u - ua -o.io 


06+0- 04 

00+0- 01 
u,uu -0.00 

00+ ' 02 
u ' uu -0.00 

00+ 02 
u - uu -o.oo 

04+o- 04 

u,u ^-0.04 


30+0- 01 
u.ou_ 01 

44+0-01 
u, ^-0.01 

34+ 01 
u -°*-o.oi 

42+0- 01 
u '^ z -0.01 

47+0-02 

u '^' -0.02 


1036.6/776 
1328.5/1007 
825.5/776 
896.2/776 
1308.2/1009 


1.34 
1.32 
1.06 
1.15 
1.30 


MCG-6-30-15 


93/07/09 
93/07/31 
94/07/23 
97/08/03 

Q7 /OS /07 


2 m+ 01 

z - ui _0 .01 

1 QQ + 0.01 

1 qo+0.01 

i «n+ - 01 

i 77 +0.01 

'-■"-o.oi 


fiK +0.03 
°- 65 -0.03 

70+ 04 

67+ 01 
u -°'-o.oi 

oqi+0.03 
u - Mi -0.03 
70+O.O3 
u -'°-0.03 


90+O.O2 
u,zo -0.02 
49+0.03 
u '^ 3 -0.03 
12+ 001 

n7+o- 03 

u - u ' -0.03 

15+0- 03 
ulo -0.03 


oo+o- 00 

oo+o- 00 
u - uu -o.oo 

011+0' 002 

u.ui -0.0Q2 

n 009+0.004 

u.uaz_ 004 

n36+o o 03 

U.U3D_o 004 


1203.7/775 
1110.9/772 
2861.5/1013 

1972.6/1013 
i fi/iq a /1 ni 


1.55 
1.44 
2.82 
1.95 
1.00 


MR 2251-178 


94/04/02 
93/1 2/07 
93/12/14 
93/12/19 
93/12/24 


i O1+0.02 
1 ' 81 -0.02 
, -c+0.02 
i -'°-0.02 
1 67+ ' 03 

1 63+ - 04 
1 -"°-0.04 
i ,-7+0. 03 


40+ 07 

41+ 06 
U -* -0.05 

40+ ' 08 
u - w -0.09 

5+ ' 1 

u -°-o.i 

41+O.O8 


q 17 +0.05 

io+o- 04 

U ' iU -0.04 

o.o9t° :S 

0-066±S:g? 

oo+o- 06 
u - uu -o.oo 


000+ 003 
u.uuu_ 000 

ooo+o- 004 

u.uuu_ 000 

000+0' 007 
u.uuu_ 000 

oi+o- 02 
u - ul -o.oo 

OOO+O' 004 
U.UUU_ 000 


829.3/823 
821 7/867 

Oil -L . 1 / OU 1 

768.0/750 
628.1/696 
634.4/660 


1.01 
0.95 
1.02 
0.90 
0.96 


NGC3516 


94/04/02 
95/03/11 
95/03/12 


i 71 +0.01 

1 - 1 1 -0.01 
i 79 +0.03 
^•'^-0.02 
i 70+O.OI 
1 -' 8 -0.02 


n 74+0.03 
u ''^-0.03 

64+ ' 05 

u -° -0.05 

64+ ' 04 
u '°*-0.04 


39+0-02 

u - oy -0.03 

21+0- 05 
U ' Zi -0.03 

16+0- 04 
u ' 1D -0.04 


0+°'° 
0.05+ 01 
0.0281°,;°°? 


924.3/772 
1025.6/1005 
1079.1/1002 


1.20 
1.02 
1.08 


NGC4051 


93/07/27 
96/06/27 


9 07+O.OI 
z -°' -0.01 
9 910+O.OIO 
z.zia_ 0Q4 


01+ 03 
u - ul -o.oi 

00+ 01 

u,uu -0.00 


55+0- 03 
u - oo -0.03 

25+0- 01 
U-^O-0.01 


0+°'° 
+ °'° 


1320.4/748 
1746.6/1010 


1.77 
1.73 



NGC5548 93/07/27 


1 92+ ' 01 
1 - a -0.01 


n 90+0.03 

U '^ S -0.03 


°.MtS:S 


0+°'° 


829.8/772 


1.07 


96/06/27 


1 91+0.02 

1 - al -0.02 


0.363llg° 4 


0.08l°'° 3 


047+0' 005 
U - U£ " -0.005 


1156.1/1011 


1.14 


96/06/29 


1 89+°'° 2 
1 - OM -o.oi 


65+ ' 04 


000+"' 009 
u.uuu_ 000 


qqi+0-006 
U.UOl_ 005 


1249.5/999 


1.25 


96/07/01 


1 84+ ' 02 
1 - oy: -0.02 


69+°'° 5 


016+ 04 

U.U1D_ 02 


041+ ' 007 

U.U^1_ 007 


1096.4/1002 


1.09 


96/07/03 


-, 7 r-+0.03 
1 -'°-0.02 


5fi+0' 04 


01+ 05 
u ' ul -0.01 


03+ 01 


995.5/984 


1.01 


96/07/04 


1 77+O.O2 
' -0.01 


5fi+°'° 2 
U-O0_o 05 


00+°' 04 
u,uu -0.00 


03+ 01 
V-Vd-0.00 


1105.1/1008 


1.10 


98/06/15 


1 94+0-01 
^ -0.01 


n 90+O.O2 
{hZS -0.03 


00+ 02 
u - uu -o.oo 


046+°'°° 4 
U.U1D_ 004 


988.4/774 


1.28 


98/06/20 


1 944+0.004 
1 - a ^ i -0.004 


n 901+0.009 

U '^ 81 -0.015 


000+°'°° 6 
U.UUU_ 000 


n^i+o- 001 
u.uoi_ 001 


2115.6/774 


2.73 


98/07/07 


1 70+O.O2 

t- '8-0.01 


16+°'° 3 


01+°'° 4 
u ' ul -0.01 


03+ 01 


1006.4/902 


1.12 


99/01/19 


1 84+°'° 3 


24+°'° 7 
u - z -0.06 


03+ 05 
U.UO_ 03 


06+ 01 
u - UD -o.oi 


697.0/664 


1.05 



Table 2. Spectral fitting absorption results for the phenomenological model, as defined in Section 3.2 of the main text. Column 
3 shows the best fitting power law photon index. Column 4 and 5 give the optical depths at threshold of the OVII and OVIII 
absorption edges respectively. Column 6 reports the best fitting column density of intrinsic neutral absorbing material (placed at 
the redshift of the source). Column 7 gives the goodness of the fit parameter and the number of degrees of freedom (dof). Finally, 
column 8 presents the Red. x 2 °f the fit. 



photoionization code Cloudy C90.04 (Ferland et al. 1998). 
Initially, a power-law ionizing continuum has been used, of 
photon index F (determined by the phenomenological model 
above), and luminosity between 2-10 keV, L2-10 equal to 
10 43 erg s _1 , incident onto a shell at inner radius r = 10 16 
cm of density n — 5 \ 10 9 -10 13 cm -3 . Therefore, the ion- 
ization parameter: 

£2-10 



(4) 



varies from 0.01 to 20 erg cm s 1 . The column density was 



varied from N n 



10- 



-10^ 



NGC 4051 is the 



only object in our sample that has F > 2 and also the only 
source presenting absorption features at energies larger than 
1 keV . It required larger values for the ionization param- 
eter: 1-50 erg cm s _1 ; hence the range of values covered 
by the density is n = 5 x 10 10 -10 n cm ~ 3 . We will be as- 



suming a constant density cloud i n our modelli ng, which 
may not be true in every case (e.g. Binette 1998). A black 
body has been included for MCG-6-30-15, NGC 3516, NGC 
4051, and NGC 5548 as given in Table 3. Recall that spectral 
variability of MCG— 6-30-15 strongly argues for a multi-zone 
warm absorber (e.g. Otani et al . 1996, Morales, Fabian and 



Reynolds 2000, [Lee et al. 200l[ ). The simultaneous UV/X- 
ray studies of NGC3516 by Kriss et al. (1996a,b) also imply 
the existence of a complex stratified ionized absorber. The 



Ch an dra observation of NGC 4051 (Collinge et al. 2001) has 
revealed a two-zone WA for this object. Thus, one-zone pho- 
toionization models can be only regarded as a useful param- 
eterization of the spectral data. 

All the observations for each source were examined us- 
ing Cloudy models of the type described above. The results 
of these fits are shown in Table 4. The data still show devi- 
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MCG-6-30-15 



Date of 
Observation 


Photon index 

r 


kT 

(keV) 


TO 7 


TO 8 


Intrinsic Nn 
(10 22 cm- 2 ) 


X 2 /dof 


Red. x 2 


Lb (10 43 
erg s _1 ) 


Lb/L2-10 


1993/07/09 
1993/07/31 
1994/07/23 
1997/08/03 
1997/08/07 


1 96+ 02 
i on+0.02 
1 - oo -0.04 
1 m+O.Ol 
1,y -0.01 

L -°°-0.01 
i 01 +0.01 
1 ' 81 -0.02 


09+ 01 

H+ 001 

in+ 01 
u - lu -o.oo 

0.080_ 004 
10+ ' 02 


o-6il° :g5 

n 65+ 003 

n 7s +0.02 
U.78_ 002 

70+ ' 04 
u -' u -0.04 


10+ 003 
19+ - 04 
H+ 001 
04+ 01 

n 13+003 


03+ 03 
H+ 03 
05+ ' 02 

n+ 001 

-0.01 

10+ 001 
u - lu -0.03 


797.5/767 
785.8/764 
1286.3/1005 
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1994/06/07 2 -01±o.oi °- 13 ^aoi °- 22 ±ao3 - 12 ±o.02 °- 02 ±aoi 1145.4/1005 1.14 0.5 15.0 
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Table 3. Spectral fitting results for the phenomenological model including a black body component for MCG— 6-30-15, NGC 3516, 
NGC 4051 and NGC 5548. Column 2 gives the best fitting power law photon index. Column 3 presents the best fitting black body 
temperature. Column 4 and 5 give the optical depths at threshold of the OVII and OVIII absorption edges respectively. Column 6 
reports the best fitting column density of intrinsic neutral absorbing material (placed at the redshift of the source). Column 7 gives 
the goodness of the fit parameter and the number of degrees of freedom (dof). Column 8 presents the Red. x 2 of the fit. Column 9 
reports the total (bolometric) luminosity of the best fitting black body and Column 10 shows the ratio of this bolometric luminosity 
to the 2-10 keV luminosity. 



ations from this best fitting at ~ 1 keV . This presumably 
signals a breakdown in one of the assumptions of this model 
(such as the one- zone nature or the chemical abundances of 
the plasma). 

An illustration of the incident and transmitted contin- 
uum as calculated by our models and the confidence contour 
plot for two observations of MCG-6-30-15 and NGC 5548 
is illustrated in Fig. lQ. In Table 5 the values for L a b s for 
each of the 30 observations is reported. 



3.3 Results and comparison with other estimates 



3.3.1 Mean values, errors and comparison with other 
methods 

We have substituted the results of L a b s and iV wa from each 
observation in Eqn. 2. The value of the black hole mass, Mbh 
for each observation versus the 2-10 keV luminosity, L2-10 
is presented in Fig. 2 and reported in Table 6. There are dif- 
ferences in the results for Mbh obtained for each observation 
that can be accounted for within the errors, except for one 
observation of NGC 3516^J A possible cause for these dif- 
ferences could be the dynamical evolution of the WA, since 
we would be proving different material for different obser- 
vations. In particular for the case of the NGC 3516 UV and 



t The incident and transmitted continuum as calculated by our 
models for al l the observations as well as all the contour plots are 
presented at [http : //«ww-xray . ast . cam, ac ,uk/~rm/BH_mass^ 



t See Table 4. The 2-10 luminosity for this source in the 1994 
observation is approximately 2 times larger than in those observa- 
tions in 1995. Also note how our photoionization model does not 
reproduce very well this very first observation, giving a x 2 = 1.91. 
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Source 


Date of 
Observation 


log 7V wa 
cm " 


( erg cm — 2 s — 1 ) 


Intrinsic iVjj 
(10 22 cm" 2 ) 


£-2-10 
(10 43 erg s" 1 ) 


X 2 /dof 


Red. 

x 2 


IC4329a 
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Table 4. Spectral fitting results using the one-zone photoionization model. Columns 3 and 4 give the best-fitting column density 
and ionization parameter 5 of the warm photoionized plasma. Column 5 gives the best fitting column density of intrinsic neutral 
absorbing material (placed at the redshift of the source). Column 6 reports the goodness of the fit parameter and column 7 the 
Red. x 2 of the fit. 



X-ray abso rber, such evolution ha s been reported (Koratkai 
et al. llS9| and [Kriss et al. 1996b|). 



It is worth noting that all the MR2251-178 observations 
were taken within one month and for this particular object 
the consistency of the results is remarkable. There is also 
a variation in L2-10 with time for all the sources (i.e. the 
source gets brighter or fainter). These changes in the X-ray 
luminosity do not seem to relate to the values obtained for 
Mbh in a simple manner. 

For several of these AGN, there is extensive obser- 
vational evidence for a multi-zone WA. For the case of 
MCG-6-30-15. time variability of th e OVIII edge (e.g. 



Otani et al. 199t , Morales et al. 200C ) and t he detection 



of different ionization species ( Lee et al. 2001 ) support the 
multi-zone hypothesis. And very recently, for the case of 
NGC 4051, jCollinge et al. (200l| ) have resolved two dis- 



tinct absorption systems. We have not taken into account 
the multi-zone nature of the WA in the present work. 

Another element that might have some effect on our 
results is that below approximately 1 keVSISO and SIS1 
spectra since approximately late 1994 have been showing 
an increasing divergence from each other, and from the GIS 
data, towards lower energies. This is in the sense that both 
SIS0 and SIS1 efficiencies below 1 keVhave been steadily 
decreasing over time. This loss in efficiency is currently not 
corrected by any of the software and a detailed analysis of 
how this could affect our results is discussed in Section 3.3.2. 

Once the value of the black hole mass was calculated 
for each observation, we have obtained the mean value for 
each source. For that, those observations with \ 2 > 2 have 
been excluded. Our result for NGC 4051 might not be sig- 
nificant giving the high value of \ 2 f° r both observations 
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MCG-6-30-15 NGC5548 
2 1994_07_23 2 1996_06_27 




Energy (keV) Energy (keV) 



MCG-6-30-15 




21.6 21.61 21.62 21.63 21.64 21.3 21.35 21.4 21.45 

log N H (cm -2 ) log N H (cm -2 ) 

Figure 1. Top panels: example of the incident continuum (top line) and transmitted continuum (bottom line) as 
predicted by our models for the 1994/07/23 MCG-6-30-15 and the 1996/06/27 NGC 5548 observations. The area 
enclosed by the two lines gives the value for the flux absorbed by the WA clouds. Bottom panels: contour plots for 
the two observations in the plane ionization parameter, £, against column density, JV W a at the 68, 90 and 99 per 
cent confident level for the two parameters of interest. 



(although our result is consistent with other methods esti- 
mates). The errors only take into account the uncertainties 
in Labs (calculated from the dispersion in the value of L a b s 
for each observation) and iV wa (obtained from the fits). The 
mean values for Mbh are presented in Figure 3 together with 
other independent meth ods estimates: reverberation map- 
ping ( |Wandel et al. 19991) for IC4329a, NGC 4051, and NGC 
5548; X-ray variability ( |czerny et al. 200(i| ) for MCG-6-30- 
15, N GC 4051, and NGC 5 548; and velocity dispersion re- 
sults ( |Gebhardt et al. 2000| ) for NGC 3516 and NGC 4051. 
Our estimates for all sources cluster around 10 7 Mq. Despite 
the simplicity of our model, our results reproduce within an 
order of magnitude other independent estimates. 



3.3.2 The case of NGC4051 

At low energies, the X-ray continuum of NGC 4051 has been 
observed to be domina t ed by a variable soft e xcess (e.g 



in strength, pre cisely the rever se of the behavior seen in 
MCG-6-30-15 ( [Reynolds 199^ ). Thi s object has very re - 



Turner fc P ounds 1989; Guainazzi et al. 1996|) . Guainazzi 



et al. (199( ) also argued for significant variability of the 
OVII edge while the OVIII edge remained roughly constant 



cently been observed with Chandra ( Collinge et al. 2001) 
and they clearly detect a two-zone WA. 



There are only two observations for NGC 4051 and both 
of them, when fitted to our photoionization models, give a 
X 2 ~ 2. This indicates that our model breaks down for this 
particular object. In Section 3.2 it was already noted that 
NGC 4051 was the only object in our sample that has V > 2 
and also the only source presenting absorption features at 
energies larger than 1 keV . 

We have investigated if the reason for this disagreement 
could be the divergence between SIS0 and SIS1 spectra and 
from the GIS data, since approximately late 1994. To test 
this hypothesis we have analysed each detector individu- 
ally. Table 7 reports the results and Fig. 4 gives a graphical 
representation of them. Surprisingly the only result that is 
inconsistent with the rest is one given by GIS2 in the first ob- 
servation. Apart from this single value, the other estimates 
are consistent with those obtained for the case of the full 
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Figure 2. Black hole mass, Mbh, in solar units for the 30 observations of the 6 AGN (the names of the objects 
appear at the top left corner of each figure) versus the 2-10 keV luminosity, L2_io- The solid line corresponds to 
L = Z/Eddi the dashed line to L = O.lZ/Eddi an d the dot-dashed line to L = O.OlZ/Eddi where L = 101/2-10- 



analysis, favoring the singular nature of NGC 4051 as the 
cause for not fitting a single ionization parameter model. 



4 THE INFLUENCE OF THE UNDERLYING 
CONTINUUM ON OUR ESTIMATE 

The emitted spectrum of Type 1 galaxies in the extreme 
UV/very soft X-ray band is poorly known although the pho- 
tons are highly effective at ionizing warm absorbing gas. 
Study of the effect of the excess soft X-ray emission on our 
results is therefore of great importance. The object we have 
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Figure 3. Our results for the mean black hole mass (stars) of the six sources versus the mean 2 to 10 keV luminosity, 
Z/2— 10- Other estimates (reverberation mapping (Wandel, Peterson & Malkan 1999), marked as squares for IC4329a, 
NGC 4051, and NGC 5548, X-ray variability (Czerny et al. 1999) , marked as triangles for MCG-6-30-15, NGC 
4051, and NGC 5548, and velocity dispersion (Gebhardt et al. 2000), marked as circles for NGC 3516 and NGC 4051) 
have been plotted at our mean L2— io- Reminder: our results are upper limits on the black hole mass and 
are not marked as such to avoid confusion. The solid line corresponds to L = Z*Eddi the dashed line to L = O.lLEddj 
and the dot-dashed line to L = O.OlLEddi where L = IOL2— 10- 



NGC4051 




2 




Figure 4. Values obtained for the BH mass of NGC 4051 when 
the 4 different ASCA instrument are independently fitted for the 
1993 and 1994 observations of our sample. The dot-dashed line 
corresponds to L = O.OlLEddi where L = IOL2— 10 • 



selected to study the influence of a different underlying con- 
tinuum has been NGC 5548. The reason for that has been 
the availability of a model for its intrinsic spectrum (i.e. af- 



ter the effect of Galactic absorption at the soft excess and 
reddening a t the optical/UV spectr um has been removed), 
as given by Magdziarz et al. (199c). Also for this particu- 



lar object there are 10 different observation in our sample, 
which allow us a better comparison with our results. 

Fig. 5, left panel, presents a montage of the incident 
and transmitted continuum as calculated by our photoion- 
ization models for 4 observations of this object. The different 
soft excess for each observation generated a different level of 
ionization in the WA and therefore the spectrum for differ- 
ent observations present different absorption edges (e.g. for 
observation 1998/07/07, labeled 9 in Fig. 5, the absorption 
edges for Hi and He n at at 0.0136 keVand 0.0544 keV, 
respectively, are very prominent). However, since the soft 
X-ray spectrum has to be reproduced by the models after 
the ionizing radiation has passed through the WA, the black 
hole mass estimate through the ratio Labs/N-m would reg- 
ulate itself to changes in either quantity (i.e. higher values 
for the luminosity absorbed demand higher values for the 
column density and a strong EUV excess means that He 
is more highly ionized). As seen in Table 6, the errors can 
account for the differences in the results. 

The comparison between the BH mass estimates ob- 
tained using our und erlying continuum, {M^n /Mg))cont.a., 
and those using the Magdziarz ct al. (1998|) spectrum, 
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Figure 5. Left panel: Incident (dashed line) and transmitted (solid line) continuum for the 1 (1993/07/27), 6 
(1996/07/04), 8 (1998/06/20), and 9 (1998/07/07) observations of NGC 5548, as calculated by our photoionization 
models. This figure shows the differences in the incident continuum for these observations, as well as graphically 
presenting the different amounts of luminosity absorbed by the WA. The absorption edges for Hi and Hen at at 
0.0136 keVand 0.0544 keV , respectively, become very prominent for some observations. Right panel: Comparison 
of the results for the BH masses using our underlying continuum, (AfBH/A^©)cont.a, versus those obtained using 
Magdziarz et al. (1998) continuum, (M w /Mq) 

cont.b- The dotted line marks the equality of both values. The 1 
(1993/07/27), 6 (1996/07/04), 8 (1998/06/20), and 9 (1998/07/07) observations are labeled. 



(AfBH/Afo)cont.b, is illustrated in Fig. 5, right panel. The 
former results appear to be lower than the latter, but still 
the differences are accounted for within the errors. 



5 ANALYSIS OF FURTHER TERMS IN THE 
EQUATION OF MOTION 

5.1 The drag term 

This section discusses the possibility of a drag force acting 
on the WA clouds, and also presents some problems that the 
cloud model raises. For the purposes of estimating the black 
hole mass, the gas that is responsible for the WA is assumed 
to exist as isolated clouds which are likely moving through 
an intercloud (IC) medium. These clouds could be confined 
by a hot IC medium. A two-phase model, i.e. cold clouds em- 
bedded in a hot IC medium, for the broad and narrow line 
regions in quasars, BL R and NLR, respectively, has been ex- 



tensi vely studied (e.g. 



Krolik et al. 1981 



Mathews & Murray 



1987). The major objection to the model is the cloud confine- 



ment problem (if the clouds are not continuously formed, an 
external mechanism is necessary) . Another problem is cloud 
break-up by drag forces. The difficulties with the model, 
especially in relation to the hot medi um temperature, are 



discussed in 



Fabian et al. (1986) and Mathews & Ferland 



(1987). The latter reference discusses th e dynamical impl 



cations of the two-phase model. Also in Elitzur & Ferland 



(1986) the stability of these clouds is addressed. Our WA 



cloud model will be subject to all these caveats. The IC 
medium could have a very significant effect on the dynamics 
of the clouds through its drag and the resulting hydrody- 
namic instabilities. Whether the WA clouds can move rela- 
tive to a IC medium without suffering gravitational, tidal, 
shear, or Ray leigh- Taylor instabilities is beyond the scope 
of this work. We assume that magnetic fields internal to the 
clouds must be important. 



First we note that ignoring the drag force justifies hving 
chosen our estimates as upper limits. If the density of this 
external gas is large enough then the WA would be subject to 
a drag force that would cause us to overestimate the strength 
of gravity in the dynamics of the WA. The upper-limit of the 
black hole mass is found by assuming that the radiation force 
on the WA, T^ad, is larger than the gravitational restoring 
force, -Fgrav Therefore, the drag force will be important if it 
is comparable to F rac i: 

^t-^abs 



Pov A. 



(•») 



The density of the IC medium at which drag will be impor- 
tant can then be related to the velocity of the WA through 
the medium (which we take to be the observed radial veloc- 
ity): 

Labs ,„\ 

no = - — 5 5-, (6) 

where p is the mean mass of particles in the intercloud 
medium (assumed to be a pure hydrogen and helium mix- 
ture). Making use of the average Labs for each object, we 
plot in Figure 6 this critical density as a function of the out- 
flowing velocity of the WA. A distance of r = 10 16 cm was 
assumed, and the range in v was chosen to cover the range 
that is typically observed in WAs. If the density of the IC 
medium is below the line then drag is not important; alter- 
natively, drag will be important if the IC density lies above 
the line. For typical values of the observed WA blueshift, 
an IC density > 10 7 cm ~ 3 seems to be required for drag to 
be important in any of the objects. The most direct limits 
of the gas come from opacity arguments (i.e. the absence 
of soft X-ray absorption). We shall note at this stage that 
a medium with a density of 10 7 cm ~ 3 and a pathlength of 
10 16 cm has a Thomson depth of nearly 0.1. 

Another constraint on the properties of the IC medium 
is that its column density has to be smaller than that of the 
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Source 


Date of 


log L a b s 




Observation 


(erg s - 1 ) 


IC4329a 


93/08/15 


42.6 ± 0.3 




97/08/07 


42.3 ±0.3 




97/08/10 


42.4 ± 0.3 




97/08/12 


42.2 ± 0.3 




97/08/15 


43.0 ±0.3 


MCG-6-30-15 


93/07/09 


42.60 ± 0.05 




93/07/31 


42.70 ± 0.05 




94/07/23 


42.59 ±0.05 




97/08/03 


42.67 ±0.05 




97/08/07 


42.64 ± 0.05 


MR 2251-178 


94/04/02 


42.42 ± 0.06 




93/12/07 


42.37 ± 0.06 




93/12/14 


42.33 ± 0.06 




93/12/19 


42.37 ± 0.06 




93/12/24 


42.24 ± 0.06 


NGC3516 


94/04/02 


42.608 ± 0.001 




95/03/11 


42.609 ± 0.001 




95/03/12 


42.606 ± 0.001 


NGC4051 


93/04/25* 


42.85 ± 0.05 




y4/UD/ U i 


a o on _i_ n r\x. 


NGC5548 


93/07/27 


42.3 ±0.2 




96/06/27 


42.4 ±0.2 




96/06/29 


42.5 ±0.2 




96/07/01 


42.2 ±0.2 




96/07/03 


42.6 ±0.2 




96/07/04 


42.5 ±0.2 




98/06/15 


42.3 ±0.2 




98/06/20* 


42.4 ±0.2 




98/07/07 


42.5 ±0.2 




99/01/19 


42.1 ±0.2 



Table 5. Luminosity absorbed by the WA. Column 1 
gives the name of the source, column 2 indicates the date 
(year/month/day) of the observation and column 3 the values of 
the luminosity absorbed, L a b s - Those observations marked with 
an asterisk, *, have a value of x 2 > 2 when the data is fitted to 
our photoionization model (see Table 4). 

WA, A^wa; otherwise the IC medium could be identified with 
the WA itself. 



Source 


Date of 
Observation 


logL 2 _io 

(erg s- 1 ) 


logM BH /M Q 


IC4329a 


93/08/15 
97/08/07 
97/08/10 
97/08/12 
97/08/15 


43.87 ±0.07 

44.01 ± 0.07 

43.88 ± 0.07 
43.93 ± 0.07 

44.02 ± 0.07 


7 O + 0.4 
'• Z -0.4 
7 q+0.4 

7 9+0.4 
'• z -0.4 
7 !+0.4 
'• x -0.4 
7 q+0.3 

'■^-o.s 


MCG-6-30-15 


93/07/09 
93/07/31 
94/07/23 

n"7 /no /no 

97/ US/ 03 
97/08/07 


43.10 ± 0.06 
43.03 ± 0.06 
43.09 ± 0.06 

a o nn I n r\c 

42.99 ± O.Uo 
42.96 ± 0.06 


6 90+ 07 
R s1 +0.09 

7 nr+0.06 
7 OK+0.07 

'• U0 -0.06 


JVLH 2251-178 


c\ a Id a /no 

94/U4/02 
93/12/07 
93/12/14 
93/12/19 
93/12/24 


a A m I n n/" 

44.91 ± U.Oo 
44.96 ± 0.06 
44.90 ± 0.06 
44.88 ± 0.06 
44.80 ± 0.06 


6 - 9 -o.i 
6 9 +0 1 

6 9+ ' 1 
70 +o.i 
'• u -o.i 


NGC3516 


94/04/02 
95/03/11 
95/03/12 


43.4 ±0.1 
43.2 ± 0.1 
43.2 ± 0.1 


6.88tS;°i 

'■ UO -0.02 

7 04+ a02 

'• u *-0.02 


NGC4051 


93/04/25* 
94/06/07 


41.53 ± 0.07 
41.62 ± 0.07 


6 47+O.O8 
6 45+ 06 


NGC5548 


93/07/27 
96/06/27 
96/06/29 
96/07/01 
96/07/03 
96/07/04 
98/06/15 
98/06/20* 
98/07/07 
99/01/19 


43.69 ± 0.07 
43.87 ±0.07 
43.83 ± 0.07 
43.75 ± 0.07 
43.68 ± 0.07 
43.77 ±0.07 
43.85 ± 0.07 
43.91 ± 0.07 

43.70 ±0.07 
43.74 ± 0.07 


6 9+ 02 

O.3_ 2 

7 1+0.2 
'•^-0.2 
7 9+0.2 
' -0.2 
7 9+0.2 
'• Z -0.2 
7 9+0.2 
'• z -0.2 
7 9+0.2 
'• z -0.2 
7 9+O.2 
'■ Z -0.2 
7 q+0.2 
' -°-0.2 

7 4+0-2 

'•^-0.2 
7 1+0.2 
'• x -0.2 



Table 6. Black hole mass, Mbh, in solar units for the 30 ob- 
servations of our study. Column 1 gives the name of the source, 
column 2 indicates the date (year/month/day) of the observation 
and columns 3 and 4 give the 2 to 10 keV luminosity, L2— 10 > 
and the black hole mass, Mbh, m solar units. Those observations 
marked with an asterisk, *, have a value of x 2 > 2 when the data 
is fitted to our photoionization model (see Table 4). 



5.2 The acceleration term 

The acceleration term has to be positive for our estimates 
to stand up as upper limits for the BH mass, and not lower 
limits as it would be for the case of negative acceleration, 
or just measurements of the BH mass for the case of zero 
acceleration. The dynamics of WA clouds in pressure e qui- 
librium with external magnetic pressure are studied by 



Ch 



elouche & Netzer (2001). Their principal conclusion is that 



such clouds can be accelerated to high velocities by means 
of radiation pressure acceleration. The velocity profiles for 
two generic cases of clouds (shell-like clouds with constant 
mass and constant density clouds) are presented in their fig- 
ure 10. For most external pressure profiles, the acceleration 



results are positive for typical values of WA velocities. Still 
for some cases a negative acceleration can be possible. 

Finally we would like to mention that for the WA 
temperature (around around 10 5 K ) , the plasma would be 
highly supersonic. Therefore, pressure forces would not be 
relevant for the velocities considered (e.g. several hundred 
kilometers per second). 



6 DISCUSSION 

In this paper we present a new method for obtaining an up- 
per limit on an AGN BH mass by means of the comparison 
between the radiative and gravitational forces acting on its 
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Date of 
Observation 


ASCA 
apparatus 


log L 2 -10 
(erg s - 1 ) 


log M BH /M Q 


93/04/25 


SISO 
SISl 


41.6 ±0.1 
41.3 ±0.1 


6 4 +0 ' 4 
"■^-0.4 

D '°-0.4 




GIS2 


41.6 ±0.1 


'■°-0.5 




GIS3 


41.7 ±0.1 


fi C S+ 11 
°-°-0.4 


94/06/07 


SISO 
SISl 


41.6 ±0.1 

41.7 ±0.1 


fi 5+ ' 4 
D '°-0.3 

fi 

°-°-0.3 




GIS2 


41.6 ±0.1 


64 +0.5 




GIS3 


41.6 ±0.1 


6 4+ - 4 



Table 7. Absorbed luminosity and BH mass for NGC 4051. First 
column reports the year/month/day of the observation, second 
column gives the ASCA when the 4 ASCA instrument used in the 
fit. Finally columns 3 and 4 give the 2—10 keV luminosity, L2—10, 
and the black hole mass, Mbh, in solar units. 
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Figure 6. The density at which any external medium becomes 
dynamically important is plotted as a function of the radial veloc- 
ity of the WA. For typical values of the velocity, a density > 10 7 
cm ~ 3 is required for drag to have an effect. A distance of 10 16 
cm — 3 for the WA has been assumed. 



cloud, Labs, and the column density of the WA, A'wa. Both 
of them are obtained when the soft X-ray data is fitted to 
our single zone photoionization models for the WA. Despite 
the simplicity of our method, our estimates for the upper 
limit of the BH mass fall within one order of magnitude of 
other independent methods results (see Fig. 3). 

The only requirement of our technique is the presence 
of outflowing (i.e. blueshifted) material which kinematics 
is dominated by radiation pressure. UV absorbers could 
be used to determine the AGN BH mass using the same 
method. Narr ow UV absorption feat ures are shown in many 
Type 1 AGN ( Crenshaw et al. 1999 ) and in most cases they 
appear blueshifted with respect to the systemic velocity. 
Their typical velocities are ~ 1000 km s _1 . 

One uncertainty in our method is the presence of other 
non-radial forces that affect the dynamic of the WA, and 
therefore change its radial equation of motion (e.g. centrifu- 
gal or magnetic forces) . Further studies of the WA dynamics 
could possibly provide more information on which forces are 
acting and therefore constrain more tightly the determina- 
tion of the BH mass. Also if the mass is well-determined by 
some other means, application of our method can lead to 
information on the drag and acceleration forces acting on 
the WA. 

Metallicities different than solar would have an effect on 
our estimates. Given that Mbh oc N~ a oc Z~l , where Z wa is 
the metallicity of the warm absorber, then solar and super- 
solar abundances maintain our estimates as upper limits. It 
would be only in the case of subsolar abundances when our 
upper limit could be lower than the real value of the black 
hole mass. Note again that if the black hole mass can be 
accurately determined, our method could be used to study 
the metallicity of the warm absorber. 

It is also worth noting that in the calculation of the ab- 
sorbed luminosity, L a bs, we have used the fact that the WA 
is optically thin material. This holds for the ASCA data used 
in this paper (more considerations on the determination of 
L a bs would be required if we were to work with higher reso- 
lution data). Also better quality data would lead to a more 
accurate determination on the errors of our upper limits. 

One of the main goals of AGN studies is to measure 
the BH mass and a myriad of techniques has been explored 
over the last two decades. The use of WAs as vehicles to 
determine the AGN BH mass reinforces the SMBH paradigm 
and opens up a powerful and simple technique for measuring 
AGN BH masses. 



WA. For that, we have adopted a cloud model for the WA in 
which the radiation pressure dominates the dynamics of such 
clouds. No further assumptions regarding the geometry, lo- 
cation and nature of the WA are necessary. This method has 
been applied to five Seyfert 1 galaxies (IC 4329a, MCG— 6- 
30-15, NGC 3516, NGC 4051 and NGC 5548) and one radio- 
quiet quasar (MR 2251-178), which cover a range of almost 4 
orders of magnitude in X-ray luminosity. The calculation of 
the upper limit of the BH mass required the determination 
of only two quantities: the luminosity absorbed by the WA 
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